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Abstract. The evolution of the comoving cosmic merger- 
rate density of neutron star binaries nc{z) is calculated 
using a distribution of their merging times provided by 
population-synthesis computations of binary stars. We 
adopt an exponential law for the star formation rate 
with various timescales for different morphological types 
of galaxies. For elliptical galaxies also an initial burst of 
star formation, lasting one Gyr, is considered. The result- 
ing nc{z) of most models agree with the form nc{z) oc 



For cosmologies with no vacuum energy, the bright- 
ness distribution of burst intensities expected for a uni- 
form source population is consistent with the BATSE dis- 
tribution if the limiting rcdshift is about unity (e.g. Der- 
mer 1992; Mao & Paczyhski 1992). However, cosmological 
time-dilation effects in the BATSE sample indicate that 
the dimmest sources should be located at z ~ 2 (Nor- 
ris et al. 1995). If this is the actual limiting redshift, a 
source population with a comoving rate density of the 



(1 + ^) 



1.5-2 



for 



2, which has been proposed for the form nc{z) oc (1 -I- z)'^ is compatible with the BATSE dis- 



source population of 7-ray bursts. Assuming a standard 
candle luminosity, the computed brightness distribution 
is consistent with the BATSE results if bursts at the peak 
flux threshold, P = 0.4 photons cm^^ s~^, are located at 
a limiting redshift of 1.9 to 3.3. Progenitors of the sys- 
tems producing 7-ray bursts at small redshift (bright) are 
likely to host in spiral galaxies and star forming regions 
whereas these at high redshift (dim) reside mainly in el- 
liptical galaxies. The location of a burst may be up to 
~ 1 Mpc away from the host galaxy. 

Key words: stars: binaries - stars: evolution - stars: neu- 
tron - cosmology: theory - gamma rays: bursts - gamma 
rays: theory 



1. Introduction 

The recent discoveries of optical counterparts of the 7-ray 
bursts GRB 970228 (Groot et al. 1997) and GRB 970508 
(Bond 1997) and the measured emission-line redshift of 
z = 0.853 (Metzger et al. 1997) for the latter provide ev- 
idence for their cosmological origin. A possible source for 
7-ray bursts (GRBs) may be the merging of two neutron 
stars in a close binary (Blinnikov et al. 1984). 



tribution for 1.5 ^ P ^ 2 (Horack et al. 1995). Therefore, 
it is worthwhile to examine whether the comoving merger 
rate density of neutron-star binaries evolves in a similar 
fashion. 

By means of population synthesis computations for bi- 
nary stars the merger rate of neutron-star binaries [here- 
after (ns, ns)] can be computed (see Portegies Zwart & 
Yungclson 1998 and references therein). Using this ap- 
proach Lipunov et al. (1995) computed the evolution of 
(ns, ns) mergers as a function of redshift and log N — log S 
distributions of GRBs for a cosmic population which con- 
tains galaxies with a constant star-formation rate and also 
galaxies with initial bursts of star formation in different 
proportions. Totani (1997) computed the evolution of the 
GRBs rate density from (ns, ns) mergers in a model based 
on the observationally determined history of cosmic star 
formation and in a model derived from detailed galaxy 
evolution. Totani assumed a distribution of merging times 
f{tc) oc t^^ for (ns, ns) systemad. We go a step further 
by adopting for all galaxies exponentially decreasing star 
formation rates (SFR) with different timescales depend- 
ing on galaxy morphology (e. g. Sandage 1986). We apply 
the distribution of tc from model B of Portegies Zwart & 
Yungelson (1998, henceforth PZY98) which provides the 



^ Sahu et al. (1997) made similar computations for "cosmic" 
star formation history and a fixed merger time of 3 x 10^ yr. 
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best fit to the expected birthrate and orbital parameters 
of the Galactic population of high- mass binary pulsars. 

2. Model 

The annual (ns, ns) birthrate is proportional to the SFR 
in the galaxy \E'(i, r), where r is the timescale for star 
formation. The SFR adopted is a decreasing exponential, 
which appears when the SFR is supposed to be propor- 
tional to the gas content without taking into account the 
gas ejected by stars (e.g. Bruzual & Chariot 1993): 

*(igal> ^) = exp(-tgal/T), (1) 

where the age of the galaxy is given by igai ^ t — tp 
and tp the time at which the galaxy was formed. At any 
time t the (ns, ns) merger rate /i(t, t) has a contribution 
from systems that are formed at different epochs t — tc in 
the history of that galaxy. It can thus be expressed as a 
convolution integral of the birthrate of stars \E'(t, r) and 
the distribution of (ns, ns) merging times f{tc)- 

l-t-tF 

fi{t, r) = /xo / f{tc) *(i - tc, t) dtc, (2) 
Jo 

The results of the population synthesis computations 
give f{tc) and the normalization coefficient /ip. The func- 
tion f{tc) from PZY98 is approximated by a Gaussian 
with X = log(tc/^gai(0)) as a parameter, the maximum at 
Xo — —2 and with a=l (tgai(O) is the age of the galaxy at 
z — 0). We assume that the distribution function f{tc), 
normalized to unity, is time-independent and the same 
for all galaxies. For the normalization of /i(<, t) we re- 
quire that in a reference galaxy similar to the Milky Way 
(type Sb, Mj-ci = 2 x IO^-'^Mq, with a current astration 
rate of 4M0 yr~^) the current rate of (ns, ns) mergers is 
2 X 10^^ yr^^, which is the merger rate obtained for the 
Galaxy by PZY98. 

For simplicity, we split the Hubble sequence into three 
types of galaxies: E to SO, Sa to Sb and Sc to Sd. For our 
selected mixture of Hubble types, two sets of characteris- 
tic star formation timescales are used (cf. Table 1). This 
parameterization is rather simplistic; the fraction of E/SO 
galaxies may differ considerably from 20% (Dressier 1980) 
and the star formation history in spirals of the same mor- 
phological type may be a function of their mass (Gallagher 
et al. 1984). If indeed most stars in the Universe formed 
in dwarf star-burst galaxies (e. g. Babul & Ferguson 1996) 
at 2: '--^ 1, this may affect our results considerably. 

Star formation is assumed to occur continuously ac- 
cording to Eq. in all galaxies. We also investigate the 
case of an initial burst of star formation during the first 
Gyr in E-SO galaxies and no star formation thereafter. 
The latter models are denoted as sets lb and 2b. 

The comoving rate density can be related to /i(z,T) 
via the B-band luminosity density from the Universe Lb = 



Table 1. Adopted parameters for galaxies of different mor- 
phological types: timescale of star formation r;; contribution 
to the B-band luminosity of the Universe Ci (Phinney 1991); 
mass to blue-light ratio M/ Lb (after Lipunov et al. 1995 and 
Guiderdoni & Rocca-Volmerange 1987). 



Type 


n (Gyr) 


Ci 


{M/LB)^ 




set 1 set 2 




[h X {M/Lb)q] 


E - SO 


1 1 


20 % 


10 


Sa - Sb 


4 6 


40 % 


5 


Sc - Sd 


10 15 


40 % 
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1.9 X 10** h Lqb Mpc~3 (Efstathiou et al. 1988) where h = 
i?o/100 kms-^Mpc^^: 

"=W = ]^ E c.(|^) ^M-,r.). (3) 

Here the summation is taken over all morphological types 
and Ci is the contribution of each type of galaxy to the B- 
band luminosity of the Universe, M/Lg is in solar units 
(see Table 1). 

The comoving rate density depends on the cosmolog- 
ical parameters (/i, flo, Ao) and on the redshift of for- 
mation zp, via the relation for the age of galaxies rel- 
ative to redshift: igai(z) = t{z) — t{zp). We consider 
two cosmological models with zp = 5 and no vacuum 
energy: {h, fto) = (0.5, 1.0) and (0.75, 0.2) for which 
^gai(0)=12.16 and 9.92 Gyr, respectively. 

3. Application 

Figure 1 shows the redshift dependence of nc{z) for sev- 
eral of our models. We normalize the merger rate density 
to unity at z=0. The gray shaded area demonstrates that 
the results are basically consistent with the estimates de- 
rived from the GRB source distribution as observed by 
BATSE (Horack et al. 1995); the computed merger rate 
falls roughly in the allowed region for z 2. Qualita- 
tively, the behavior of ndz) is easily understood from 
Fig. 2 which shows the relative contribution of galaxies 
of different types to the (ns, ns) merger rate (for set 1). 
At high z it is dominated by early-type galaxies with ini- 
tially a high time-averaged SFR. In set lb (upper solid 
line) the sudden drop at z ~ 2.5 reflects the cessation of 
star formation in ellipticals in one Gyr after formation and 
the relatively short typical (ns, ns) merger time of ;$ 1 
Gyr (PZY98). The more gradual decline of ndz) in set 
1 reflects the continuous star formation in both elliptical 
and spiral galaxies. Comparison of sets 1 with 2 suggests 
that the evolution of ndz) is better described by the mod- 
els with longer star-formation time scales. The lower solid 
line uses the same merger-time function f{tc) as Totani 
(1997). This model predicts a lower (ns, ns) merger rate. 
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Fig. 1. Evolution of tlic comoving (ns, ns) merger rate den- 
sity as a function of rcdsiiift, normalized to unity at z=0. The 
shaded area corresponds to the allowed region according to 
nc{z) oc (1 -f z)^ with f3 in the range 1.5 to 2. The upper solid 
(dotted) line is for set lb (set 1) with fgai(0) = 12.16 Gyr. The 
dashed lines give the rates for set 2 with fgai(0) = 12.16 Gyr 
(upper curve) and 9.92 Gyr (lower curve). The lower solid line 
has to be compared with the dotted line and corresponds to 
the merger time distribution f{tc) oc with a lower cut-off 
at 0.02 Gyr (like in Totani 1997). 




Fig. 2. Relative contribution of galaxies of different types to 
the (ns, ns) merger rate density. The current galactic age is 
igai(0)=12.16 Gyr and set 1 of the SF timescales is used. 



The spectral form of the burst similar to that observed is 

adopted: 4>(S) oc E''^ cicp{-E/ Eq), where £;o=350 keV 
is a characteristic energy. For comparison, we use the ob- 
served integral brightness distribution from the BATSE 
3B catalog in the energy range 50-300 keV measured at a 
timespan of 1024 ms (Meegan et al. 1996). Figure 3 shows 
the expected brightness distributions computed for sets 
of parameters 1 and lb, superimposed on the BATSE 
data. The curves are normalized at the peak flux threshold 
P=0.4 photons cm~^ s~^. Figure 4 provides the results of 
a Kolmogorov-Smirnov test of the BATSE 3B catalog to 
the results of our computations. Only data above the peak 
flux of 0.4 photons cm~^ s~^ are used to avoid threshold 
effects. The highest confidence level (CL) is obtained for 
sets 1 and lb if the limiting redshift zqa at the peak flux 
threshold is ~ 2.4 and 3.0 respectively. For other values of 
the star formation timescales and cosmological parameters 
zoA is found to range from 1.9 to 2.7 (with a CL > 80%) 
for the models with a burst of star formation in elliptical 
galaxies (models lb and 2b). Without an initial burst of 
star formation zqa ranges from 2.9 to 3.3. Similar results 
are found for Eq in the range 300 - 400 keV. 

Note that for higher limiting redshift, Fig. 4 shows the 
existence of other possible fits with a lower CL (i.e. model 
lb, Zo,4 = 3.8). The primary peak at zqa = 2.4 corre- 
sponds to the first change of the slope of the comoving 
rate density (Fig.l). For higher redshifts, the sudden in- 
crease of the merger rate would require the same behavior 
of the BATSE data for consistency. As a consequence, the 
secondary peak at zqa = 3.8, although providing a good 
fit for peak flux values near the threshold, tends to de- 
part more and more with the data for higher values of 
P (Fig. 3). Therefore, the redshift range related to an as- 
sumed initial burst of star formation in elliptical galaxies is 
likely to be beyond the actual limiting redshift for GRBs. 

Finally, it has to be pointed out that the models can 
hardly reproduce the peak flux values of BATSE's faintest 
bursts (sec also Totani 1997). As shown by Reichard & 
Meszaros (1997), this feature results from the assumption 
that GRBs are standard candles. 



The upper and lower dashed lines in Fig. 1 (set 2) 
demonstrate the effect of the age of galaxies on the co- 
moving merger rate, i. e.; the influence of the cosmological 
parameters. 

Figure 2 gives the relative contribution to the (ns, ns) 
merger-rate density for each selected subclass of galaxies. 
The majority of events which are potentially detectable 
by gravitational- wave observatories (LIGO/VIRGO) are 
located in early- type spiral galaxies. If GRBs originate 
from (ns, ns) coalescence the dimmest bursts are expected 
to be hosted in elliptical galaxies. 

Following the standard procedure (e.g. Horack et al. 
1996) we compute the number of bursts N{> P) with a 
peak flux greater than P. We assume that bursts are stan- 
dard candles and the intrinsic luminosity doesn't evolve. 



4. Conclusion 

The computed relative merger rate as a function of red- 
shift is in agreement with the detected rate of 7-ray bursts 
up-to the limiting redshift of BATSE. For the models 
where elliptical galaxies experience an initial burst of star 
formation the occurrence rate of (ns, ns) coalescence in- 
creases suddenly by more than a factor two at a redshift 
of 2.5, which is beyond our current detection limit. The 
synthetic log N — log P distribution is compatible with the 
observations down to the completeness limit of 7-ray cata- 
logues. The same is true for models based on observation- 
ally inferred cosmic star formation history (see e. g. Sahu 
et al. 1997) and for models based on galactic evolution 
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Fig. 3. Cumulative counts of bursts N{> P) of the BATSE 3B 
catalog (histogram) and model curves computed for parame- 
ters shown in the figure (for tgai{0) — 12.16 Gyr) 
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Fig. 4. Redshift zoa necessary to achieve agreement with the 
BATSE 3B log - log P distribution. The solid line corre- 
sponds to set 1 and the dotted line is for set lb. In both cases, 
the current galactic age used is tgai{0) = 12.16 Gyr. 

(Totani 1997). They, however, predict a different behavior 
if star bursts occur at high z. 

The absolute value of the (ns, ns) merger rate is found 
to be ~ 100 times larger than the GRBs frequency. Es- 
cape from this conundrum is obtained if the opening angle 
of the observed phenomenon is a few degrees, which is con- 
sistent with fire-ball models where leptons are converted 
into bulk barionic motion (Meszaros & Rees 1992). 

If 7-ray bursts indeed originate from (ns, ns) mergers 
and our model for star formation is correct, interesting 
implications follow. The progenitors of bright GRBs or 
bursts at low redshift (z ^ 1.5) most likely belonged to 
early-type spiral galaxies whereas the progenitors of the 
dimmest bursts and those at high redshift (z ^ 2) were 
located in elliptical galaxies. The majority of the bursts 
of gravitational waves in this model are expected to orig- 
inate from early- type spiral galaxies. If most stars in the 
Universe formed in dwarf star-burst galaxies, a substan- 
tial fraction of the parental population may originate from 
them. Suggested by current observational data the decline 



in the star formation rate beyond z « 1.2 (e. g. Connolly 
et al. 1997), if real, will show- up as a turn over in the 
logiV — logP distribution for GRBs. 

As noticed by Tutukov & Yungelson (1994) and con- 
firmed by PZY98 (their Figs. 6 and 8), space velocities 
of (ns, ns) binaries may well exceed escape velocities of 
dwarf as well as giant galaxies and they may travel up to 
~ 1 Mpc before coalescence. Thus, a significant fraction 
of the sites of GRBs may not be directly associated with 
star forming regions. 
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